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ABSTRACT 


The equations for estimating the performance of an axial flow 
turbine are examined to determine the efficiency and other parameters 
when very high power outputs per stage are required. It is indicated 
that turbine stages should be designed for peak efficiency or for 


high power output, but that compromise designs are not advisable. 


A survey of the immediately available literature on losses in 
turbine stages at the U. S. Naval Postgraduate School is outlined. 
An attempt to find and correlate specific loss data with theoretical 
relations was not successful. Data is fragmentary, theoretical re- 
lations are few, and the correlation of the available data on profile 
losses, tip clearance losses, effect of trailing edge thickness and 


secondary losses is poor. 


An axial flow single stage turbine was designed for transonic 
flows and a digital computer program was written to evaluate the 
off-design performance of this turbine. The particular relations for 
supersonic expansion after a blade row are included in the program 
as well as the effects of entry shocks into the rotor. Initial 
indications are that the program should be useful for evaluating 


turbines designed by the procedure shown. 


The computer programs included are written in FORTRAN language 


for the CDC 1604 Digital Computer. 
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THE PERFORMANCE OF HIGH POWER OUTPUT AXIAL 


FLOW TURBINES UTILIZING FLOWS IN THE TRANSONIC REGIME 


I. Introduction 


The increasing need in space applications for turbines of high 
power output from a minimum physical volume and weight requires more 
knowledge of efficiency, losses and performance estimation of these 


turbines. 


If the required power output and space limitations are predominant 
then efficiency must suffer. As a means of Showing the degree of 
efficiency sacrificed, the equations which determine the power output 
and efficiency of an arbitrary turbine were programmed for digital 
computer solution in terms of the head coefficient. The head coef- 
ficient for peak efficiency was determined and then the head coef- 
ficient was increased in increments and various parameters affecting 
the turbine performance were computed at each increment until the 
efficiency dropped to an arbitrary cutoff point. The data are plotted 
in non-dimensional form and simple equations are found to represent 


the data. 


Of paramount importance in the design of a turbine is the esti- 
métion of the losses due to boundary layer formation, mixing effects, 
clearances and secondary flows. This area is explored in the hope 
of finding a more accurate and precise method of presenting these 


losses. A survey of the immediately available literature is included. 





A particular single-stage turbine is designed for supersonic 
flow leaving the stator blade row and a relative subsonic flow as 
seen by the rotor. The equations for determining the effects of 
after-expansion behind a blade row and entry shocks entering a 
blade row are programmed for the computer in addition to the 
regular performance parameters. The rpm and pressure ratio across 
the turbine are then varied to show the effects of demanding extreme 
power outputs from a turbine so designed. A performance map of the 


turbine is developed. 





Il. Development of High Head Coefficient Data 


A. Basic Equations 
In Ref. 1 the equations are developed to establish the design 

parameters for turbine stages. These equations are obtained by 
assuming that the flow through the turbine is represented by the 
conditions on a mean stream surface and the primary parameter is 
taken to be the “head coefficient". The head coefficient relates 
the total isentropic drop of the stage from the inlet stagnation 
pressure to the discharge static pressure to the kinetic "energy" 


of the peripheral speed at the rotor entrance and is defined as 


Ay = oe (1) 


In order to relate the theoretical expressions to values found 





in engineering use, some expression for the losses involved must be 
used. A velocity coefficient (Y) is used as a measure of rotor 


efficiency and is defined as 


We 2h - 


This velocity coefficient has been found through experience to depend 
primarily on the deflection of the flow in the rotor (Af). Vavra 
has established a mean curve through data from a variety of sources 


to be 
2.28 4A.97 
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and this equation is used to represent the losses for purposes of 


evaluating the efficiency at each step of the calculation. 


Following the procedure outlined in Ref. 1, arbitrary values 
were chosen for the angle of the flow leaving the stator COr 
reheat factor (f), and the carryover coefficient of kinetic energy 
into the next stage (D,): the radius ratio (R,/R,); the meridional 
velocity ratio (V_,/V_,), the velocity coefficient in the stator (Y), 
and the carryover coefficient of kinetic energy from the stator into 
the rotor (Pp): Then the degree of reaction Ge) is selected and 
the head coefficient is increased in increments. The degree of 
reaction is the fraction of the isentropic enthalpy change through 
the turbine stage which occurs in the rotor. At each increment of 
head coefficient the particular values of velocity coefficient (VY), 
internal efficiency (7): leaving loss coefficient (k,) > ratio of 
meridional absolute velocity to peripheral speed (Vy /Up) > relative 


velocity ratio (Wo /W,)> and the relative and absolute angles (OX; 


Px 2 are computed. 


A similar computer program using the same equations was written 
to find the peak efficiency for a given entering angle and degree of 
reaction, and then all the other parameters are computed. The values 
are also found for higher head coefficients such that the efficiency 


is 0.5% less than optimum and then 1.0% less than optimum. 


rr ay, ee 7 


oe 6 ee 4h ee tm i ome e 





B. Computer Program 


The equations to be solved by iteration until A ae! are; 


anna: Pine 
(1 ‘Paced, -A® lh (4) 
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Once this condition ( vn — WY ) is satisfied the values of the 








other parameters are found from 


T= ie ae oe AEE) 
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The equations were coded in Fortran language, Ref. 2, and solved 
on the Control Data Corporation 1604 computer at the U. S. Naval 
Postgraduate School. A copy of the program is presented in Appendix III. 


The resulting data printouts are included in Appendix V as Table I. 


The program finding the peak efficiencies, efficiencies 0.5% less 
than optimum and efficiencies 1.0% less than optimum is also in 
Appendix III, with data in Appendix V as Table II. It should be 
noted here that all equations are completely coded in Fortran language 
so that any selected numerical values may be introduced in the be- 
ginning of the program. Any of these arbitrarily chosen values may 


be changed and the equations of the main program remain unchanged. 


C. Results 
Representative data from Tables I and II are presented in Fig. l. 
The regularity and progression of the data curves indicate that non- 


dimensionalizing the efficiency and head coefficient parameters might 
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he profitable. When all the efficiencies for one entering angle and 


one degree of reaction are divided by the optimum efficiency for that 


set of conditions (71 q 


ope? and the head coefficients are divided 


by the optimum head coefficient (k, /k; . a the curves of Figs. 2, . 


3, 4 and 5 result. 


A closer inspection of these plots reveals that for each degree 
of reaction, one curve suffices for all four entering angles. 


The cubic equation 


Kae ; ae i i 
M7 ) = -.000223 |=] + .008345{ ~S——_|_ -. 12068 { $8] 41. 195 
opt keys 
is opt is opt 1S opt 


* 
represents all the data for r = 0.0 for all entering angles from 60° 


to 75°. 


* 
Similarly as an approximation to all the data at r = 0.25, the 


equation is 


a7 Kis : is e is 
C717 yo.) ~-000223 | #2 | 4.00770 {48——|" -.1098 | 8 — | 41. 
opt k k, k, 
is opt is opt 1s opt 
* 
and for r = 0.50 
Ks 3 Ki. 2 Keg 
W7 ) = -.000125 | ———— +.00554 | ——~-———_- -.09414 | —————]_ +1.097 
opt Ke k, k. 
1S opt 1s opt 1S opt 


{ne maximum error in the range of entering angles presented here is 
approximately 8 parts in 130. This is probably accurate enough for 


preliminary engineering estimates based on the equations in Ref. 1 








used in the computer programs. As will be seen later, the loss data 


is not known to this accuracy. 


The curves show that, roughly speaking, increasing the head 
coefficient to five times the optimum causes a reduction of 30% of 
the obtainable peak efficiency. Increasing the head coefficient to 
ten times the optimum causes a reduction of 45% of the obtainable peak 
efficiency. This indicates that if it is necessary to go to very high 
head coefficients to obtain a certain required work output, the rate 
of change of efficiency with increasing head coefficient is greatest 
from optimum to five times optimum head coefficient. The rate of 
change of efficiency decreases as the head coefficient is further 


increased. 


Table I of Appendix V can be used to establish the other stage 


parameters after the head coefficient range of interest is chosen. 


Figs. 1 through 5 are not the off-design performance of a 
particular turbine, but they are the locus of points of an infinite 
number of different turbine designs. Each point on each curve 
represents a turbine designed for that condition. Therefore, if the 
requirement is that the turbine be designed for a high head coef- 
ficient, the loss in efficiency need not be the absolute controlling 
factor, and the design may as well be based on extreme high values 
of head coefficient. These investigations show that it is advisable 
to either design the turbine for peak efficiency of maximum output 
but do not compromise both by operating at only two or three times 


the optimum head coefficient. 





III. Loss Investigation 


Since the design of a turbine depends on the loss assumptions 
made, the use of an average velocity coefficient based on general 
considerations of flow deflection in the rotor is not a very precise 
method of estimating the losses. This was the method used in the 
preceding section of this thesis. A better way should be found and 


the following considerations are basic to the problem. 


A. Loss Components 

It seems logical to class the losses encountered in turbines 
into four particular categories; 

1) Profile loss -- the friction loss associated with the 
rermation of the boundary layer on the blades of the cascade. 

2) Tip clearance loss -- the pressure and energy loss associated 
with the leakage flow between the tips of the blades and the casing 
ot the rotor, including the effects of trailing vortices at the hub 
where there is zero clearance. 

3) Trailing edge loss -- the loss due to mixing effects behind 
the blades of the rowas caused by the finite thickness of the blade 
trailing edge. 

4) Secondary loss -- the loss due to the vortices which are set 
up in the flow in the curved channel between blades giving a velocity 


component normal to the desired direction. 





In the axial flow turbine, the closely packed blades and large 
flow deflections usually require the flow through the turbine to be 
considered from the point of view of "channel flow' between blades 
instead of using airfoil theories. As will be seen in the subsequent 
discussion, this secondary flow causes the eae ae problem pre- 


venting the accurate assessment of the losses. 


Looking at the flow between two blades of an arbitrary cascade, 
Fig. 6, the velocity distribution is as shown for the mid-span of 
the blade channel and it can be seen that a static pressure gradient 
is developed between surfaces 0-1-2 and 0'-1'-2'. The centrifugal 
forces on the fluid being turned cause an increase in the static 


pressure near the 0'-1'-2' wall and decrease the static pressure 





near the 0-1-2 wall. 


Fig. 6 Flow Between Blades of an Arbitrary Cascade 
=|0= 





Assuming non-viscous flow outside the boundary layer requires the 

flow to increase velocity near the 0-1-2 convex wall to maintain the 
total pressure a constant (P = p+ 4 pw). Conversely, the flow 
velocity must decrease near the 0'-1'-2' concave wall outside the 
boundary layer. In the boundary layer at the top and bottom of the 
channel, the loss of velocity (and total pressure) due to friction 
causes the flow direction to be dominated by the pressure gradient 
imposed (assuming the static pressure is constant through the boundary 
layer). Therefore, as the flow progresses from station 0-0' to 
station 1-1' to station 2-2', the fluid in the upper and lower 


boundary layers tends to flow from 0' toward 1 and from 1' toward 2. 


This change of direction of the flow in the boundary layers 
results in the formation of two vortices superimposed on the main 
fiow leaving the channel at station 2-2'. In Fig. 7 the flow patterns 
in the upper and lower boundary layers are shown as well as the 


induced vortices rotation direction at the exit plane. 


This vorticity results in velocity components developed perpen- 
dicular to the desired exit velocity. Since a velocity component 
normal to the desired direction cannot be recovered usually, most of 
the useful energy in these components is lost in friction. This is 
the so-called "secondary loss''. For a rotating row of blades, the 
blade tip clearance also has a marked effect on the secondary flow. 
The leakage from the high pressure to the low pressure side of a blade 


encourages an additional trailing vortex to form a reaction blading. 


-ll- 





For impulse blading the flow over the tips which is not deflected also 
enhances vortex formation. -For cascades of infinite blade height 
(two-dimensional only) there is no secondary flow since the mechanism 


depends’ principally on the annulus boundary layers. 


inte cm oe 





Fig. 7 Flow in the Boundary Layer of a Cascade Channel 


With the physical parameters of blade height, trailing edge 


thickness, tip clearance, profile shape, and spacing given, it should 


be possible to make a reasonably exact estimate of losses. 


B. Survey of the Literature 
The immediately available literature at the U. S. Naval Post- 


graduate School was surveyed to find the particular loss breakdown 
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used by the various sources, the theory and representations offered, 
and any precise loss data published. The survey is limited in scope 
and critical since the particular desire is to locate specific loss 
theories and data. The results of this survey are presented here 

in the order of reference books, general texts, and then periodicals 
and reports. The NACA/NASA reports provide a wealth of data on the 
efficiency, or total loss coefficient, of specific turbine designs 
from tests. Very little information was found for loss components, 


however, in the frame of reference considered here. 


1. In Ref. 1 an equation is derived for the determination of 
a flow coefficient in a rotor (Y) as a function of the change in 
flow angle of the relative flow. This equation is a fitted curve 
through data from several sources and a variety of blade shapes and 
operating conditions and is the same equation used in the preceding 
section of this thesis. In Fig. 8 is shown the general shape of the 


curve and the derived equation (3). 


A loss coefficient can be defined as 


Ym /—y* 


Considering the case of 48 = 0, we see the loss coefficient 


(11) 


¥ =A - y ? represents those losses in the cascade which are 
present when the flow is not deflected. As the flow is deflected 
an amount AB, the reduction in YW indicates an additional loss. 
The additional loss will be taken as the secondary loss. The 
secondary loss coefficient is then 


no wo * 


= 





Curves shown in subsequent discussions of secondary losses will now 
Show a loss coefficient curve based on this equation plotted as a 
dashed line and labeled "VAVRA". This curve will act as a general 
reference and will show "average" conditions so a comparison of each 


method can be made. 
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Secondary Loss 





_ 228 4,97 
= 09 = ——4f~— 
: oat (80-48 


As 


Fig. 8 Flow Coefficient W as a Function 
of Turning Angle in a Cascade - 


2. A veritable fountainhead of theory, design and ideas are 


presented in Ref. 3. This basic work by Stodola touches on almost 


all the subjects considered important now. The important works of 


wee 





others up until 1927 are included to show their measurements and de- 
velopments. The losses through fixed rows of blades are given in 

the form of the so-called velocity coefficient Y . the ratio of actual 
velocity to theoretical velocity. The measurements of the Institution 
of Mechanical Engineers of England and others all show good correlation 
for flows up to sonic speed for stators with an axial entering angle 
and an exiting angle of 70°-78° from the axial. For supersonic 

speeds the data show some scatter due to difference in measuring 
apparatus and some basic differences in design of the convergent- 


divergent passage. 


Secondary flows in bends are explained but there is wide vari- 
ation in the available data for friction losses in bends. No pub- 
lished data were available for the losses in rotating blade rows. 

The data and experiments on fixed rows of impulse rotor blades were 
of limited application and only show general trends. 

For design purposes, velocity coefficients and loss coefficients 
are given in general terms only, i.e., "velocity coefficient for 
stators = 0.975-0.92 for long and short nozzles respectively". The 


curve presented by Stodola for the velocity coefficient for rotating. 
blade rows as a function of turning angle is included in the data 


used for the determination of the equation (3) by Vavra. 


The tip clearance losses had been explored in a limited fashion. 
These losses were accounted for by an empirical formula for the de- | 


crease in internal work as a function of clearance and blade height, 
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¥ ¥ ona en (13) 
where k = radial clearance and h is the blade height corresponding 
to zero clearance. This relation, credited to Anderhub, is 


plotted on Fig. 9 and noted as (A) : 


The steam turbine design procedures presented are strongly 
dependent on making the new design a small modification of an already 
designed unit which is in service and for which the efficiency 
is known. The peak efficiencies were recognized to be in the regions 
of k from 1.0 to 5.0, as shown in the previous section of this 


is 


paper. 


3. The axial turbine section of Ref. 4 gives a lucid and brief 
explanation of the flows in such machines but the loss data and 
secondary flow considerations are based completely on Ref. 22, 23 


and 24 which are considered later in this paper. 


4. Since books written to be used as school texts of necessity 
try to cover the complete spectrum of fluid mechanics, theory, and 
design, only the most cursory reference is made to specific losses 
in turbine stages. Refs. 5, 6 and 7 are representative of books of 
this type. Ref. 5 defines a work recovery coefficient as the ratio 
of useful work extracted by the rotor to the total kinetic energy 
of the flow leaving the stator, times the stator flow efficiency. 
This work-recovery coefficient variation with the velocity ratio 
(U/V,) is shown for impulse and reaction conditions and no other 


explanation of the nature of the losses is offered. The charts 


Bibs 





presented are based on data from Ref. 16, in which report the losses 
in a blade row are presented only as a function of stagger angle and 


incidence. 


In Ref. 6 the author makes the statement that no adequate 
criterion is available for the difference in the flow conditions 
between an efficient and an inefficient hydrodynamic machine of a 
given type. Any losses and efficiencies of turbines offered are 
tied to airfoil theory or definitions of losses by Stodola and 
others. Secondary flows and radial flow of the boundary layer are 
mentioned only from Ref. 17 where the meaSurements were made on 


widely spaced blades on a rotating blade row at low speed. 


A large list of references on the work in gas turbines is 
given in Ref. 7 for the data available at that time. Principally 
the work in Ref. 16 below is mentioned in connection with turbines 
but only general terms are offered and no specific loss derivations 


are offered. 


5. The most comprehensive effort to make a breakdown of losses 
and show specific data found in a textbook is that of Ref. 8. This 
book takes selected data from Refs. 15, 16 and 23, discussed later 
in this paper, and presents plots of the losses due to leakage, 
profile losses, losses from turbulence and wall friction (including 
secondary loss effects) and incidence losses. The discussion is 


brief, however, and no new data or theory are presented. 
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6. An attempt to mathematically tackle the problem of the 
secondary flows in bends is presented in Ref. 9. The equations of 
motion of an ideal fluid are used and frictional effects are ignored. 
An induced drag pressure drop coefficient is defined to account for 
the power loss due to secondary flow. The method is applied to the 
flow in the corner of a 3' x 4' wind tunnel which has two turning 
vanes to simulate a blade row and the method shows the proper trends. 
A quantitative value derived for the loss is computed and found to 
be about twice the measured loss. No other example of the theory 


being applied is available. 


7. A qualitative investigation of the flow in boundary layers 
and wakes of blade rows is offered in Ref. 10. An excellent discussion 
of the general nature of these flows is accompanied by the statement 


that a quantitative theory for these effects is not available. 


8. An analytical method of estimating turbine performance 
described in Ref. 11 is based on losses due to incidence angle and 
losses due to flow of the fluid through the blade row. A blading 
loss parameter is used to represent the losses but only a general 
range of values of the parameter are given. No method is presented 
for determining the magnitude of the blading loss parameter before 


the turbine is built. 


9. The general state of the art to 1948 is described in Ref. 12. 
The paper is not intended to be used as a basis for other than general 


design considerations, however, and some empirical relations presented 
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were attempts to cover both axial turbines and compressors. This 
necessarily limited the worth of such relations for large flow 
deflection turbine loss estimations. Those relations presented for 
strictly turbine use are repeated in Refs. 24 and 25 which will be 


covered in more detail later in this thesis. 


10. The losses in stators are, of course, easier to measure 
chan the losses in rotors. Some excellent tests and measurements 
are presented in Refs. 13 and 14. These show that, in general, a 
properly designed nozzle for a turbine will have a velocity coef- 
ficient (ratio of actual velocity to isentropic velocity) value of .96. 
This single value is found to be uniform over a considerable range 
of designs and conditions. Off-design values do fluctuate, partic- 
ularly in the case of supersonic nozzles but the reasons do not seem 
obscure and loss estimations for stator rows should be on relatively 


firm ground. 


ll. Ref. 15 is mentioned as a reference work in many textbooks 
on turbomachinery. This compilation of the then current design 
practices only shows the losses as a function of axial length of the 
blades of the rotor, turning angle in the rotor and incidence angle. 
The curves for loss as a function of turning angle are included in 
the data used for the determination of the equation (3) by Vavra. 


Specific detailed loss breakdowns or data are not available. 


12. The systematic investigations into secondary loss presented 


in Refs. 18 and 19 are of real interest although the test data is too 
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limited to permit any universal application. An attempt is made to 
define the two-dimensional loss by theory, measure the end loss at 
the blade root and the clearance loss at the tip, and subtract these 
from the total pressure loss measured to arrive at the secondary 
loss. The tests are performed on a particular airfoil shape in 
three solidity arrangements and three flow deflections. The inac- 
curacies and errors discussed should help any further investigator 


in this subject. 


13. A comparable breakdown of the losses into components as 
assumed here is offered in Ref. 20. This work by Markov also includes 
the experimental results of some other Russians in the field. With- 
out attempting to completely describe the development of equations 
and data presented, the pertinent data for the present consideration 


is extracted. 


The profile loss, or two-dimensional loss which would be caused 
by blades of infinite length is shown on Fig. 10 along with other 
data to be described later. The basis of the Ref. 20 data is a whole 
series of experimental measurements of widely used cascades for 
turbines in Russia. For stators, the profile loss coefficient is 


almost constant at .03 over a wide range of conditions. 


The tip clearance loss is shown on Fig. 9 as represented by two 
plots in Ref. 20. One plot is shown as the loss for a reaction stage 
and one plot is based on experimental data for the efficiency of stages 


with a small degree of reaction. The magnitudes of the degree of 
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reaction are not specified but the trend is correct. Stages with a 
high degree of reaction have a higher loss across the tip, as would 
be expected due to the pressure difference across the rotor from 


entrance to exit. 


The trailing edge loss for cases where the pressure gradient 
is small along the surface of the trailing edges of the blades is 


computed by 


» | Ae 
Ng os (14) 


Where such a pressure gradient is present the use of boundary layer 
thickness parameters is required. The trailing edge loss coefficient 


according to eq. (14) is shown in Fig. 20. 


The energy lost due to the secondary flow is treated as a 
function of the mass flow rate in the cascade. The efficiency of a 
cascade of blades with finite length is defined as 


7 - 7 | - Energy lost due to secondary flow 
(2-dim) Energy available flowing through cascade 


(15) 


The loss coefficient is defined as 


=) = (16) 


y= Fhe 7 ; Ya-asm) es = 


and so 


Energy available 


:- Y o-dim | 7 (2-dim) | ¥. | 


2 1« 





The secondary flow loss coefficient ¥ is converted to the x 


used for comparison in this thesis by 


. i x / (2-din) 


An expression for 3, is developed that is a function of the circu- 
lation around the blades, the blade width, and incompressible flow, 


Ln fy ~ Ler h 
(4/2) 


where the coefficient K depends on the expansion of the flow in the 





Fm k 


cascade. The value of K is determined from experimental data to be 


koe Pa 


refs 


he OC, 


Thus 


- Ad fz Han by i: Fan 
y= go) Ta ia 


in Fig. 11 the function ze is shown for the two-dimensional efficiency 


of Ref. 1, 
Z 
in nw % = 09262 


Values of h/c plotted are chosen arbitrarily but Ref. 20 specifies 
that the derivation holds for blades that are not excessively short. 
There is assumed to be some two-dimensional flow in the channel be- 


tween the blades at the center section. Experimental data are shown 
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tor short blades, and these data commence at ( =] = 2.0 and go to 
smaller values. So the plot in Fig. 11 is based on eq. (17) and 
is shown for a value of (2) = 2.0, the lower limit permitted. For 
( 2] =z 1.33 and 1.25 the mean coefficients are about 10% of the rep- 
resentative values based on the Vavra equation (3), but these are 
specific data for short blades not arrived at using the equation 
above. The secondary loss as computed by Markov is low compared 
to any other available data. When complete loss coefficients are 
computed in Ref. 20 the values are comparable to any other method 
however, since the tip clearance losses, trailing edge losses, 
blade height corrections, etc., used bring the total loss to the 
proper level. Any use of the Ref. 20 data must take account of 


this fact. 


Of interest also in Ref. 20 is a relation for the discharge 


flow angle of a turbine blade row. The relation is 


CoLrx, = —_+—_—_—. 

Ys (18) 
4—(A / 602 %, ) 

and it gives more accurate results than the other widely-used 


equation 


A 
flor X, = —— 
AL (19) 


which does not take account of trailing edge thickness. 
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14. The tip clearance loss as represented by Ref. 21 is also 


shown in Fig. 9. The equation derived is 


5 


Y= ——_4+—__ (20) 
Ch/ 4) + & 
where § , = 1011 + 4.667 (k/b) 
v; = -.0422 + 2.790 (k/b) 
h # blade height 
b # blade width 
k *#® clearance 
and the results are shown for two values of (h/b). 
Also given in Ref. 21 is the equation for the minimum induced 
drag due to flow in a clearance gap at the tip of an airfoil (as 
derived by Betz) to be 
z 
D =k, hh P Wea | (21) 
where K, is a form factor depending on the clearance, blade spacing 


and relative exit angle of the flow. 


Assuming the axial component of velocity to remain constant 
and drawing an arbitrary turbine stage velocity triangle (such that 
all angles are positive), Fig. 12a, the circulation (/’) can be 


computed by integrating around the closed curve in Fig. 12b to give 


(SZ aud (Wig 2 — Was ) (22) 
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Then the lift per blade is 


L=pGIudk = ped (, (Wz - hh) (23) 


The lift and drag forces exerted on a blade can be seen in Fig. 13 
as well as the axial and peripheral components of the resultant 


force F. It can be seen that for each blade channel 


=4k (2-7) (24) 


(Wes Mar) 


wa | 
a ae (25) 


| 


The secondary flow loss will be associated with the loss from the 


pressure drop in the axial direction. 


To arrive at a loss coefficient expression for the axial force 
component, it is necessary to define the loss coefficient in terms 
of the pressure loss, as shown in Fig. 14. The loss coefficient is 


defined as 


iP he = 


Substituting 
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and for incompressible flow, 


- / a 
lp = Ay +H , (27) 


lez = Athena 


and collecting terms, 


pa 
(A-P)= Pg — 3 eae (28) 





Then from Fig. 13, 


le ~ Eten fh, + Pee f, “i 


Using the above derived expressions for L and D and equating 


4k (y-z,) = Lam, + Dts feg (30) 


eives 





Ss =) — ag | (nfs orf) an fe aA, 
+ 4K fF Laos fy an coef} 


For an arbitrary value of KF | , this secondary loss is plotted in 
Fig. 15. The value of % (5 can be adjusted to change the general 
inclination of the family of curves shown to any desired value. That 


is, a lower value of K raises the curves and vice versa. 


ah} 


-9%- 


| Name 01 me Ome ate) fay +) Prem ey 
SN SE a Ge kA ert tee ee Ste ee ® 
— ee eve 6 6 we te & (ee a) 0) 
bee OF em me om Rigs Ve eee ee 





15. In Ref. 22 an attempt is made to find a mathematical 
solution to the problem of secondary flows in particular. The 
approximate relations derived are not able to satisfy the boundary 
conditions. One conclusion drawn from the derivations is that the 
secondary flows are a direct consequence of a non-uniform approach 
velocity, and, if the lift is assumed proportional to the square 
of the local velocity, then the secondary flow is independent of 
aspect ratio. These conclusions are counter to the ideas of most 


investigators in the field. 


The tests and measurements made showed odd discrepancies and 
it was felt that rotor-stator induced effects may have been the 


cause. No relation was deduced which is directly useful here. 


16. The ideas and relations presented in Ref. 23 and 24 were 
carried through and refined and elaborated in Ref. 25. All three 
references haye the same authors. Therefore, Ref. 25 will be dis- 
cussed here as being most useful, and parts of Ref. 23 and 24 will 


will be extracted as necessary. 


By far the most detailed method presently available for deter- 
mining loss coefficients is presented in Ref. 25. The loss is 
broken into component parts similar to those assumed at the begining 
of this discussion and the quantitative data offered ne based on 


tests of turbines. 
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The loss coefficients of Ref. 25 are 


P ik 


ay R2 a Loss of total head pressure 
. Peo - Po Total pressure at blade outlet-static pressure at blade outlet 
(32) 


The loss coefficients used in this thesis are 


ler —le2 








7-2  @ 
_ / « le, Pa a 
and the necessary conversion is 
ka, — 'ez 
Y = kez —~ Fz eae 
/ +- fe, — Fer sae (34) 
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For profile losses, two figures are presented, one for nozzle 
blades and one for impulse blades. These are reproduced here as 
Fig. loa and Fig. 16b. An interpolation method is proposed for rotors 
with reaction and a method of extending the profile loss to conditions 


of other than zero incidence is also shown. 


A number of relations for the loss due to tip clearance are 


presented in Ref. 25 as 


f= 3/fh 
a (Stodola, Reaction turbine) (35) 


x = 45 (-f (Meldah1) (36) 
aff a 


¥ = 26(4 
4 (50% Reaction turbine) (37) 


These relations are shown in Fig. 9. 
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Expressions are also derived for tip clearance loss based on Cy 
and C. if the ratio spacing/blade height is small and the gas turning 
angle is small. For our uses here, the gas turning angles are large. 


The loss relation arrived at is : 
2 
a ap 
Vy - 15 (2. Cod a, [coer Man X, Tan %, (38) 


where B # 0.5 for blades with clearance and B # 0.25 for blades 


with shrouds. 


Equation (38) is used to find the range of loss coefficients 
for the impulee conditions of cos X 5 = 1.00 and o, = a, - Shown 


in Fig. 9 is the band of values corresponding to x, = Xf, = 60° to 75°. 


The data available to Ainley and Mathieson led them to the con- 
clusion that the effect of area ratio is more important to secondary 


flows than is the magnitude of the gas deflection. The area ratio 


used is 
Ay (annulus area at outlet from blade row) (cos x.) 
Ay 7 (annulus area at inlet to blade row) (cos @ y) (7 


where x, = gas outlet angle measured from the axial and A: = blade 
inlet angle measured from the axial. The hub ratio (inner diameter/ 
outer diameter) was also felt to exert a strong influence. This led 
to the plot of (A,/a,)°/L1 + (ID.OD)] versus a parameter A, used in 


defining the drag coefficient due to the secondary flow as 


Ce A a [a/e) a 
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This plot is reproduced here as Fig. 17. Ultimately then, a loss 
coefficient for secondary flows is presented as 
az ' F 
KO = 4A | cone, [00d || Zink Lend 


Converting these values of i to Ke. a plot is presented as Fig. 18 


for a wide span of inlet and outlet angles. 


For trailing edge thickness losses, the assumption is made that 
the previously derived losses are for a blade of trailing edge thick- 
ness which is 2% of the blade pitch. For other trailing edge thick- 
nesses a correction is applied according to a plot which is reproduced 


heré as) Fig. 19. 


C. Results of Loss Investigation 

The results of the loss investigation will be discussed according 
to the original breakdown of loss components mentioned at the beginning 
of this section of this thesis. The specific turbine design used in 
the last section of this report will be used as the object of a 
quantitative comparison of the loss relations shown here. 

1. The profile loss estimation method presented by Ainley in 
Ref. 25 is straightforward and complete and should give good data 
for the “conventional blade shapes from which the data is derived. 
An anomaly presents itself in Fig. 10 however, when a comparison is 
made between the profile loss coefficient data of Ainley and the data 


of Markov. The blade shapes from which the Markov data are derived 
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are very similar to those of Fig. 25. It would appear that such blade 
shapes would have higher profile losses than those based on airfoil 
designs which usually are used for lower flow deflections. The data 


comparison shows this is not the case. 


ww 


Of special interest in the Markov data in Fig. 10 is the fact 
that an optimum solidity of (s/c) @ .625 to .650 is evident. Higher | 
solidities apparently block the flow thereby increasing the profile 
loss, and lower solidities encounter separation effects since the 
blades act more like individual airfoils. The low point in the curve 
corresponds to guiding the flow in a channel, retarding separation 


onset, but not blocking the flow. 


2. The tip clearance losses shown in Fig. 9 seem to give a 
large range of values depending on which curve is chosen. The sit- 
uation is perhaps more reasonable when the ratio of radial clearances 
to blade height is restricted to realistic values less than .0Q3. 
The loss estimate from the Ainley equation is obviously outside the 
flow deflection range intended for that relation. The other curves 
for tip clearance loss for reaction blading could be based on dif- 
ferent degrees of reaction corresponding to the fan aspect of the 
curves, i.e., high reaction blading would use the highest curve. 
Some reasonable estimate could be made. The only relation for 
impulse blading available is that of Markov. The Meldahl relation 
is the only one which accounts for the trailing vortex inducing a 
loss even though the tip clearance is zero, but the dependence on 


blade aspect ratio appears unreasonably strong for this data as a 
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whole. For (h/b) = 1.00, the loss at (k/h) = .02 is greater than the 


highest ordinate shown. 


3. The Markov data is all that is available for making an 
estimate of trailing edge thickness loss directly. The plots in 


Fig. 20 indicate the trends which would be expected. 


The Ainley method of modifying the complete loss coefficient 
depending on the relative trailing edge thickness seems an unfor- 
tunate way to present the data. The system is workable but indirect. 
It would be more useful for comparison and design purposes to have 


the trailing edge thickness data presented directly. 


4. The situation regarding secondary flow loss coefficients is 
the most confused of all. The general equation by Vavra is not very 
precise but certainly indicates the proper trend if the model pre- 
sented in Fig. 7 is correct. None of the other presently available 
methods investigated will even show the proper trends for high flow 


deflections. 


The Ainley relation plotted in Fig. 18 seems to have depended 
too heavily on stator data with an axial entering angle and low flow 
deflections in the rotor. Otherwise, the curves indicate that for 
a given entering angle, increasing the flow deflection decreases the 


loss. 


In Fig. 15 the relation derived from Meldahl and Betz indicates 


the proper trend for low flow deflections only. Taking the plot for 
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A = he; as an example, it can be seen that for flow deflections 
greater than Ag = 90° the indication is that increasing the flow 


deflection decreases the loss. 


The Markov relation plotted in Fig. 11 suffers the same deficiency 
pointed out for the Ainley and Meldahl curves plus abnormally low 
secondary loss coefficients throughout for reasonable blade aspect 
ratios. The data is of no value unless used specifically with Markov 
relations for all losses. In Ref. 20, when sample loss estimations 
are carried out, the total loss coefficients compare well with the 
measured values. Apparently the profile, tip clearance and trailing 
edge loss coefficients with the blade height corrections used in the 


complete method counter the low secondary loss estimates. 


5. The relations now available can be used to make a quantitative 
comparison of the component losses. The turbine designed in the 
following section will be used as a model to show the comparisons. 

The physical dimensions of the stator and rotor are as listed in 
Table IV and these are used in the relations and plots so far developed. 
The numerical results are shown in Table V for the stator and Table VI 


for the rotor. 
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Table IV 
TURBINE PHYSICAL DIMENSIONS 


ntering angle 

xiting angle 

lade height, h 

railing edge thickness, t. 
lade chord, c 

lade width, b 

pacing, s 

lade thickness, t 


ip clearance assumed, (k/h) 





Table V 
STATOR LOSS COEFFICIENT COMPARISON 

























Tip 
Clearance 
Loss 


Trailing 
Edge 
Loss 











Profile 
Loss 


Secondary 
Loss 


Author 


Javra, Ref.1 
Ainley, Ref.23 
Ainley, Ref.25 
Markov, Ref.20 
Meldahl, Ref.21 
Stodola, Ref.3 


=aLe 





Table VI 
ROTOR LOSS COEFFICIENT COMPARISON 


Tip Trailing 
Author Profile Clearance Edge Secondary 
Loss Loss Loss Loss 
Vavra, Ref.1 


Ainley, Ref.23 


Ainley, Ref.25 
Markov, Ref.20 
Meldahl, Ref.21 
Stodola, Ref.3 





* 
reaction 


As can be seen, the component losses presented in Tables V and VI 


can be combined selectively to obtain almost any desired answer. 


The loss investigation is disappointing in that no clear-cut 
answers are provided. The complete Ainley method, or the'complete 
Markov method, may produce a reasonable estimate but the correlation 


is poor at best for component losses. 


Any desire to find a way of basing a new design on component 
loss considerations is presently thwarted. A need is evident for 
an orderly testing program to find the magnitude and interaction 


of component losses. 
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IV. Design of Turbine 


As a further investigation into the characteristics of a turbine 
being used for maximum work output, a single stage turbine is designed 
for a slightly supersonic velocity leaving the stator and a relative 
Mach number of 0.8 as seen by the rotor. Then the overall pressure 
ratio is increased and the off-design performance evaluated. The 
number of calculations required are prohibitive for hand calculations 
so the equations necessary for the performance evaluation are also 


programmed in Fortran language for the CDC 1604 Digital Computer. 


A. Loss Coefficient Data 

A set of loss coefficients for the design were provided by 
Vavra from some unpublished data. These loss coefficients are shown 
in Figs. 21, 22 and 23 and will be used for the design and the per- 


formance evaluations. 


The abscissa of the plot of stator loss coefficients, Mi 5? is 
the Mach number for an isentropic expansion from the total pressure 
at the inlet to the static pressure at the discharge of the blade 
row. This can be seen in the T-s diagram of the expansion through 
a turbine shown in Fig. 24 to be the isentropic velocity - divided 


h 


by the acoustic velocity based on T When reduced, the result is 


1 is’ 


(42) 
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To properly determine the flow area through the stator, the 
loss coefficient in Fig. 20 noted "for flow rates" and shown in 
Fig. 24 as - is used. This loss coefficient is intended to 
represent the losses encountered by the flow in the channel between 
blades, but is not intended to account for the mixing losses and 
separation losses encountered by the flow leaving the blade row. 
These additional losses are included in the higher loss coefficients 
which are applied to velocity determinations. These loss coef- 
ficients are noted ''for velocities" in Fig. 21 and as ¥ in Fig. 24. 
Similar considerations are used for the loss coefficients through the 


rotor. 


The loss coefficients are corrected for blade height according 


to the curves given in Fig. 23. 


B. Flow function 
For steady adiabatic flow the stagnation enthalpy remains con- 


stant along a streamline. For stators, the stagnation enthalpy is 


za 
H = i = #, ~~ of 
(43) 
and in Fig. 24, to a different scale, H would be represented by 
To and hy would be represented by T,: A similar situation pertains 
to a relative stagnation enthalpy defined as 
ws Wes 
He = £ = 4, os (44) 


for the condition such that the peripheral velocity at the entrance 
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of the rotor equal to the peripheral velocity at the exit of the 
rotor. As shown in Appendix I, the mass flow rate for a given set 


of inlet conditions can be expressed by 








%& mel 
> 45 
yer! #7 (2) VR = 
where n is the polytropic exponent, Py is the total pressure ahead 
of the stator blade row and Py is the static pressure behind the 
blade row. 
As shown in Fig. 24, a relative total temperature and relative 
total pressure can be defined as 
Wi, 
ia / 
le, = 1, + — 
L9I&p (46) 
e i 
onl 
= kl 
IR1 = py (“24 
i, (47) 
The mass flow rate through the rotor is then 
(45a) 
Rearranging the above expressions for flow rate in non-dimensional 
terms defines a flow function, 
Dp . want .. 
48 
BAL 3 * 
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It will be necessary to investigate whether this flow function can 
be satisfied for each set of conditions imposed across a blade row. 
The choked flow condition corresponds to the maximum value of P/p 
which the blade row can accomodate. The development of the critical 
values of the flow function and pressure ratio can also be found in 


Appendix I to be 


— 
D Ze 
met] eee) (49) 


sie 
(4) = oe 
TV m+/ (50) 


A computer program to do these computations is explained in section E 





below. 


C. Geometry and Velocity Triangles 
The final design of the stator and rotor blade rows is as shown 
in Fig. 25 for the mean diameter. The velocity triangles of the 


average flow velocities at each station are as shown in Fig. 26. 


The sign convention chosen is that positive angles and positive 
peripheral velocity components point in the same direction as U, the 


rotor peripheral velocity. 
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It is convenient to define an "exit plane" and a "discharge 
plane" leaving each blade row as shown in Fig. 25 by the (e) and (4) 
designations. The design of the stator and rotor is such that the 
minimum cross-sectional area of the channel between blades occurs 
at the exit plane (e) - For maximum work output, the pressure drop 
across the stator is great enough to cause choked flow at the exit 
plane plus an additional expansion of the flow from the exit plane 
to the discharge plane. The conditions for such after-expansions 


are as shown in Appendix II based on the treatment in Ref. 26. 


In Ref. 1 a similar development is shown for a supersonic flow 
entering a blade row. Since the flow leaving the stator is super- 
sonic, it is possible under certain conditions for the rotor to see 
a supersonic flow also and proper account must be taken of the shock 


losses involved. 


D. Dimensionless Parameters 

As has already been mentioned, the number of variables to be 
handled is greatly reduced if the mass flow expression is non-dimen- 
Sionalized to form a flow function depending only on pressure ratio, 
specific heat ratio and the polytropic exponent of the expansion 


process. 


In Ref. 27 can be found a complete development of dimensionless 
parameters and referred values for minimizing the number of variables 


necessarily handled in presenting the data for a turbine performance 
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analysis. Those used in the present analysis are referred rpm, over- 
all pressure ratio, referred flow rate, and a power coefficient. The 


necessary coefficient forms are computed from 


Referred RPM * RPM/ 7 T, (51) 
Overall pressure ratio = Po/P (52) 
Referred flow rate = w sy UE (53) 
Power coefficient * HP/P Te (54) 


The referred flow rate and power coefficient are constant in this 
evaluation. It is useful to refer all the conditions at any point 
in the turbine to inlet conditions. For this reason, velocities are 


carried as V/ y 7 and W/ y i and pressures as p/P. and P,/P.- 


E. Design Computations 

In order to design the stator for choked flow at the exit plane, 
the flow function must be known for the imposed conditions. A com- 
puter program was written to compute the flow function for a given 
specific heat ratio, % , pressure ratio, P/p, and loss coefficient, ie 
As a by-product the program also computes the polytropic loss coef- 
ficient, yo This information is presented in Appendix V as Table III, 
Flow Function and Polytropic Loss Coefficients, for % from 1.25 to 
1.40, 5 from 0.0 to 0.25 and P/p from 1.02 to 6.00. The particular 
flow function data to be applied here for % = 1.37 are also included 


as Fig. 27. The computer program is included in Appendix III. 


ie 





The following data were assumed given for the design: 
Mass flow rate = 100. lbm/sec 
Total pressure at inlet = 1000 psia 


Total temperature at inlet #2 1460°R 


Specific heat ratio = 1.37 

Gas constant, R = 421.5 ft 1bf/1bm OR 
Mean diameter of stator = 23.0 in 

Mean diameter of rotor = 23.0 in 

Design rpm = 13000 


The stator blade design was provided by Vavra from unpublished 
data and is as shown in Fig. 25. It is now necessary to select the 
blade height and number of blades so the choked flow condition at 


the exit is satisfied. 


A blade height of 1.16 inches is taken as a first estimate. Then, 
from Fig. 21 the loss coefficient for flow areas at M, = 1.00 is 
found to be - = .053. From Fig. 23 the blade height correction 


is K = 0.872 and so 
f= .872 (.053) = .0462 


Using linear interpolation on Fig. 27 for this loss coefficient 


between the maximum values of D , iver, D it is found that 


cr’ 


D = .65811 


(=| = 1.86 


cr 


-42- 


i—_ ==, = a 
us & | Ff eel ee cat oe 
7 0) 0 eee eee ee - — % ot 


- te i ewe 


al wemeert weil | - Ls * fet om | 


LO + 2 et ineww & a am ee at a ~~ 


ee he «— be re “roe we a 
ry o«* - «@ 


~ Aa oy 





Recalling the flow function equation (48), it is found that 


a DQ, a — ene 


= 13.22 a 2 
_ $5811 21.0 in 


From Fig. 25 the measured minimum distance between blades at 
the exit plane is a # .448 in. and s = 1.8 in. is given. The number 


of blades is then 


d “Lr 
(35) 


In order to prevent harmonics in the wake pattern from the stator 
as seen by the rotor, a prime number of stator blades is chosen. 


In this case, z # 41 blades will suffice. For 41 blades, 


Zs 2 = 1.763 
1 


so all measured values from Fig. 25 must be scaled by the scale 


factor, 
Scale Factor #® A182 0.9794 
The blade height must now be 
A 21.0 


ot ca Gio (ete oe, oe (56) 


Recomputing p.. with a new loss coefficient based on this new blade 


height makes no measurable difference in exit area required so the 
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Stator row design is now fixed. To simplify the iteration on D a 
computer program was written to find the choked flow values of D..’ 
(P/p) and vy This computer program is included in Appendix III 


and the computed data are presented as Table VII in Appendix V. 


The flow conditions at the stator exit plane are determined 


from Fig. 24 as follows: 


fx = lL = .538 
1.86 
45 
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ly ; 
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. Ce 


Ah _ Ze _hus\(/-%.) = .1470 
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= .8458 
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When a check is now made of the Mach number at the exit plane in 


order to check for choked flow it is found that 


ye 


VEGR % 
This computed Mach number is the average Mach number of the flow at 
the minimum area. The highest Mach number at this area is equal 
to 1.00 but it is averaged with the Mach number of the flow in the 


boundary layer as accounted for by the loss coefficient. In 


2H, - 





Appendix IV the derivation of the Mach number for choked flow in a 
polytropic expansion is presented. The “critical average" Mach 


number is found to be 


/+%(#-/) G7 (Iv, 9) 


This condition will also be encountered later in the computer program 
for the flow through the rotor and proper account must be taken of 


it there. 


Using the relation from Ref. 20 for flow exit angle and the 


measured blade exit angle from Fig. 25 of Xn) = 73.9° gives 


LoL X, = lacs ils ceca = .268 
A (hy [toe Xz, ) 
(18) 
a, = 75° 


For the rotor, the blade height is permitted to be approximately 
0.30 inches greater than the stator blade height in order to permit 
the streamlines to expand and offer no resistance to flow leaving the 
stator but not have dead areas above and below the useful flow area 
of the rotor blade. The rotor blade height will be taken then as 


Ao = 1.465 in. 


In order to satisfy the requirement for the rotor to see a 
relative Mach number of 0.8, some after-expansion is required. In 
order to evaluate the proper amount of after-expansion, the equations 
in Appendix II were coded for the computer to solve for any applied 


pressure ratio between the exit and discharge planes for a given @, 
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Mach number and blade angle. This computer program is included in 
Appendix III and the solution for @ #8 1.37 and x, = 75° over a 
range of pressure ratios is presented here as Pig. 28. The computer 
program can be used for any nee rc heat ratio, blade angle, Mach 


. of 
number and pressure ratio range. 


After some iteration, it is determined that p,/P, = 0.460 should 
give approximately the correet pressure ratio across the expansion 


area from the exit plane to the gischarge plane, 


p 
1 d -460 


Now from Fig. 28 it is found that AoQ{ @0.2° and 6. % .0008. The 
total loss coefficient to be gppliedg is based on a blade height of 
1.165 inches for the stator. Tere is ¥ = »1045 from Fig. 21 and 


KR = .838 from Fig. 23 for velocity computations. Therefore 
y= (Ke) +% = .0883 


Now the velocity leaving the stator is determined with Fig. 24 


.2701 
_ (a\* -(8 = .8108 


as reference 


. 











/, 

ay a  - Jim (/-%,) ws .1892 (.9117) = .1725 
be y ~ %G 

ye ff PGR BT i 131.60 
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To find the relative angle and relative velocity of the flow, it is 


necessary to construct the velocity triangle as shown in Fig. 26. 


Gf — («ea Nie ) Den 


1G 72a. 


O = 75° - 0.2° = 74.8° 


s 34.18 
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= 92.82 
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As a check on the requirement of a relative Mach number of 0.8 as 


seen by the rotor 


ba WANE = 29887 Fico 
y \sgk7 123.9 


To design the rotor, several considerations must be kept in 


as 


mind. The flow deflection should be as great as is practicable for 
maximum work output. The trailing edge thickness will be arbitrarily 


selected as .025 inches, and the spacing and blade design will be 


chosen so the minimum cross section area of the flow passage between 
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blades will occur at the "exit" plane. The rotor blade width is 
selected as 1.375 inches and impulse conditions only will be used 


for unchoked flow, t.e., P, ™ Po: 


Recalling that a relative total temperature and relative total 


pressure are defined as 














2 
Ter hy Ww, h WY r-1) 
—— SJ — Cc — +-~ — ® 9253 
5b Fo IG OT OG tik (46) 
~ 
ee Be Cass 
RB GN GG 
(47) 
The flow function as applied to the rotor is 
Bm Vu [R 
Az fa, a (48a) 


and this expression will be used to determine the exit area of the 


rotor. For the prescribed impulse conditions 


P 
Rl | 6955 
pP, ~ 4600 ~ 171° 


and the value to be used on the abscissa of the plots for rotor 


loss coefficients in Figs. 22 and 23 


U 


i 34.18 
v, ° T3rer " - 2598 
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With this value the rotor loss coefficient for flows is found to be 
if * .819 (.192) = .1572 


Now using Fig. 27 for this pressure ratio and loss coefficient and 


using linear interpolation the flow function is 


p = 58953 


and from eq. (48a) above 


Da, * 19.13. 
=. 
Ag © 32.44 in 
An impulse rotor operating under conditions as these should have 
a ratio | 
3 
—| > 
(e] 0.7 
and 60 


s = 0.7 (1.375) = 0.962 


For the mean diameter af 23.0 inches, the number of blades on the 


rotor is then 


* 723.0) s 75.1 blades 


9 962 


Since an odd number of blades was chosen for the stator the rotor 
should have an even number of blades to reduce the possibility of 


resonance induced vibrations and the resulting metal fatigue. 


oh 





Arbitrarily choosing 76 blades 


- os = .9415 in 


and the distance between blades, a, to obtain the required area of 


32.44 square inches is 





The exiting blade angle can be seen from Fig. 25 to be 


ar Xe 
yi} 


/s2 = 71.57° 


# .3162 





oz. ™ 


The blade design itself is based on a circular arc contour of 
the concave face. The geometric relations of the design have been 
ap veened by Vavra for computer solution. Using this computer 
ea requiring solutions for 74, 76, 78, 80 and 82 blades, 
the printout of the computer solution is presented as Fig. 29. The 
values necessary to draw the blade design are shown in Fig. 30 and 


have been used to draw Fig. 25. 


The solution for 76 blades appears valid and so the other solutions 
are not used. The resulting a © 0.691. <A check on this spacing can 
be made by Brilings Rule, an early rule-of-thumb for steam turbines 


which is found in Ref. 3 to'be ; 
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{=| s 2.5 sin 2B * where s = 90 - x2 (57) 


(2) = 2.5 sin 36.86° = 1.5 


($) eae = .6667 


Also Ref. 28 presents a relation for spacing and blade angles 
based on an aeradynamic load coefficient YW, which is found to 
have a nearly constant value over a wide range of turbine designs. 


The relation is 
° ve 
Atw é wt at 
- = &, —— 4n (5A, NF) = 0.8 
Man f A- (58) 
and sa, using proper angles as defined in Ref. 28 for this case 


(2) - O.8 sin 20.407 sg agy, 


2.0 sin 18.43° sin 38.83° 


The relations by Briling and Zweifel show at least that the choice 
of 76 blades is reasonable for this design. The relation of Zweifel 
in particular is based on two-dimensional considerations and does not 
account for blade height. Some unpublished test data available to 


Vavra indicates the value of 0.691 is a good one for this size turbine. 


To complete the computations necessary to determine the power 
output of the turbine for a given set of conditions, the velocity 


triangle after the rotor is computed. 
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For no after-expansion, 4, * / 32 # Ties? , and finding loss 


coefficients in the manner previously described, 
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Using relations as before with Fig. 26 as reference, it is found 


Vee z 


= = 45.93 





4 


and the power output, with the proper conversion factors is 


ww 


HP = gd 


(U, V 


ut 7 UgVy2) (1-414) = 48,682 HP (59) 


1 


These computations are verified by the computer program. The first 


page of the printout of the computer program is presented as Fig. 31. 
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V. Computer Program 


A. General Considerations 

The computer program for the evaluation of the turbine performance 
is intended to be as general as possible within the framework of 
the turbine design method used. Any single stage turbine so designed 
can be evaluated by entering new physical dimensions, entrance con- 
ditions and loss SoStrictent data. The static pressures to be im- 
posed ahead of and behind the rotor must be programmed by the user 
arewi the main program. The loss data must be of the same format 


as Figs. 21, 22 and 23. 


B. Main Program 

1. The main program is intended to make the major computations 
and comparisons and provide control. The subroutines are provided 
to do the lengthy iterative processes and repetitive calculations. 
The flow charts for the main program and the subroutines are in- 


cluded with the program listings found in Appendix III. 


The basic concept of the computer program is as follows: 

a. With the input data given, compute the design condition 
including the after-expansion at stator discharge. 

b. Reduce the static pressure through the impulse rotor i 
increments, computing velocity triangles and power at each pressure 
increment. The shock condition for supersonic flow entering the 
rotor is computed. Also the flow function through the rotor is 


checked for the choked condition. 
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c. When the choked condition through the rotor is attained, 
the static pressure at the rotor discharge is lowered in increments 
and the after-expansion condition computed. The resulting power and 
velocity and thermodynamic conditions are computed for each pressure 
increment. 

d. The process is terminated when the loss coefficient due 
to the rotor after-expansion is reduced to zero. 

The program will also terminate if the loss coefficient due to 
after-expansion out of the stator goes to zero or if the loss coef- 
ficient due to an entry shock into the rotor goes to zero. Any con- 
tinuation of the program beyond these limiting conditions would lead 
to conditions which violate the Second Law of Thermodynamics, as 
explained in Ref. l. 

The loss coefficient reaches a maximum and then declines to zero 
or negative values as the imposed pressure ratio across the control 
area is decreased in increments for the after-expansion case. The 
exact reason for this is not known. A check was made to see if the 
axial component of the Mach number of the flow leaving the blade 
row reached M = 1.00 when the loss coefficient reduced to zero. It 
is found this is not exactly the case but that the axial component 
of the Mach number is less than one. Further work is needed on this - 


point. 


2. Some pertinent comments may be useful about the main decision 
points in the main program. The flow chart in Appendix III shows that 


the first check is for the magnitude of the static pressure between 
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the stator and rotor compared to the static pressure in the exit plane 
of the stator. For a stator designed as a converging-diverging nozzle, 
the Mach number at the exit plane may be greater than Moe = 1.00. If 
the static pressure between the stator and rotor is greater than the 
Static pressure at the stator exit plane, the shock condition is 
computed at the stator exit. If the pressures are equal, there is 

no flow deflection, and if the static pressure between the stator 

and rotor is lower than the exit plane pressure, the after-expansion 


case is computed. 


The next decision is based on the relative Mach number as seen 
by the rotor. For the supersonic case, the shock of the flow in the 
rotor entry is computed. For subsonic flow the component of relative 
velocity in the direction of the rotor blade entering angle is used 


to compute flow conditions. 


The flow function and critical pressure ratios are used to 
determine choked conditions in the rotor. As can be seen in Appendix IV 
the Mach number of the flow at the rotor exit plane is an average 
value, and the critical Mach number varies with the loss coefficient 
as shown there. Since the impulse conditions for the rotor are 
specified, some care is taken by an iteration process to find the 


P, * Po condition for choked flow before the program continues. 


Once the choked condition in the rotor is attained, the con- 
ditions forward of the rotor exit plane remain fixed. The Mach number 


of the flow at the rotor exit plane is set at Mi = 1.00 and the 
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static presqmre behind the rotor is lowered in increments. The after- 
expansion behind the rotor is then computed until the limiting con- 


dition is reached. 


C. Subroutines 
1. The subroutine called TRNGL computes the velocity triangle 
based on rpm, absolute velocity and direction or based on rpm, 


relative velocity and direction. 


2. The subroutine called AFTER computes the deflection of the 
flow due to shock or expansion of the flow out of a blade row. The 
flow angle is incremented from the blade angle and the pressure ratio 
is computed and compared to the imposed value. The process con- 
tinues until the flow deflection is found which corresponds to the 
pressure ratio imposed. The detailed equations are derived in 


Appendix II, based on Ref. 26. 


3. The subroutine for computing the flow function for a 
pressure ratio imposed across the blade row is named CPHI. The 
polytropic loss coefficient corresponding to the design flow 
loss coefficient is determined by iteration. Then the flow function 
to satisfy the imposed conditions, critical flow function, and 
critical pressure ratio are computed from the equations developed 
in Appendix I. The pressure ratio was found to be more sensitive 
in fewer significant figures than the flow function itself so 


decisions in the main program are based on pressure ratio comparisons. 


-56- 


11s oe 


>——_=—_——_«—ia —_ 


ser ‘ ot, ‘ 
= s 
- © -_= / 
2 —_ - 
—_ * 
_ ’ 


i ee ee —_ | oe 
ema i) hia 
—i—_er FF 
an os '% 1 =~ - 
. je \ah —= ‘s > 
am'¢ ° | =——" #+eenge 
oa SS ett EE = i 
o —" 7 00 —<—« Ts « > aes 
° ; tt ee .: -_' = 
| = Te mi = © 6 =e © 
- ints am i ; ttt epee 
7 = *¢ ++ 6 
: o* 4 +) Cee 6 
| -—ét— ' ° . 
* ‘ - ‘ < 7 ' 
. ' . = 
TL 







i 


—at lie © ied ogee, 





4. Subroutine BEFORE is used to compute rotor entry shocks based 
on the relations in Ref. 1 which in turn are based on considerations 
similar to those of Ref. 26 for after-expansion. Since the flow 
deflection is known, an iteration process is not required and the 
computations are straightforward. For cases which are not directly 
solvable by the basic equations, the component of the velocity in 


the direction of the blade angle is taken as the useful resultant. 


5. The CURVE subroutine is written as a means of finding a 
data point on a curve. The loss coefficient curves "for velocities" 
from Figs. 21, 22 and 23 are stored in the computer in the form of 
arrays. The subroutine writes the second-order equation of the 
appropriate three stored data points nearest to the input abscissa, 
and then finds the ordinate on the curve with the input abscissa. 
The subroutine is arranged in such a manner that the stored data 
must have equal abscissa increments and the first abscissa point 
must be zero. This means arbitrarily extending a curve such as 
Fig. 21 to M. = 0.0 whether the data at low values of M. will be 


used or not. 


D. Results of Computer Program 

1. The turbine design data given in the preceding section of 
this thesis were entered into the computer program and the resulting 
printout for the design condition only is presented as Fig. 3l. The 
pressure ratio across the turbine was then increased in increments 
(by lowering Py and Po while holding Py constant) and the resulting 


computations are included as Table VIII in Appendix V. The rpm was 
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also decreased from 13000 to 9000 in 1000 rpm increments and the 
complete pressure range imposed at each rpm. These results are 


also included in Table VIII. 


2. The non-dimensional parameters previously discussed were 
then plotted in order to make cross plots which ultimately result 


in a portion of the turbine performance map. 


3. The data immediately available from the computer printout 
permits the direct plotting of Fig. 35, Turbine Efficiency vs. 
Pressure Ratio, with the rpm as a parameter. The anticipated gradual 
loss of efficiency with increasing pressure ratio is evident. Of 
interest is the change in slope of the curve at the point where 
rotor after-expansion commences. The indication is that the increase 
in velocity due to after-expansion causes the peripheral component 
of the velocity to increase, even though a deflection of the flow 
toward the axial direction also occurs. Therefore, the rate of loss 
of efficiency is decreased, although the change is very small. A 
reduction in rpm contributes to a large loss in efficiency at any 


given pressure ratio. 


4. In Figs. 36 and 37 the most important indication is as 
follows: Fig. 37 shows that for a fixed design at a pressure ratio 
of (P /Po) = 2.174, the increase in power for an increase in overall 
pressure ratio is direct and marked up to a pressure ratio of about 
4.5 or 5.0. Further large increases in pressure ratio net very 


small gains. In Fig. 36 the corollary is shown for efficiency. The 
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initial increase in pressure ratio and power output causes a certain 
loss of efficiency but the remaining increase in pressure ratio and 


power output causes a drastic loss of efficiency. 


An interesting comparison can now be made with the high head 
coefficient data from the first section of this thesis. In that 
section, each point on each curve represents a different turbine 
design. The conclusion was that if a required high power output 
transcends the need for high efficiency, then a very high head 
coefficient should be used for the turbine design. The initial 
increase of head coefficient from the optimum value causes a large 
loss of efficiency. Once this initial efficiency loss is accepted, 
further increase in head coefficient require little more loss of 
efficiency. Now in this section concerned with the off-design per- 
formance of a particular fixed turbine design, the conclusion is: 
additional power can be obtained at a higher pressure ratio than 
design but only in the initial increase of pressure ratio. The 
extremely high loss in efficiency is encountered if the last possible 


increment of power is demanded. 


5. The turbine performance map presented as Fig. 38 serves to 
indicate the pertinent trends for the parameters investigated. The 
lines of constant efficiency form only a portion of the elliptical 
curves they become over the complete range of referred rpm and 
power coefficient. The only conditions investigated in this thesis 
are at high power coefficients. The peak efficiency for this turbine 


would occur in the lower right hand corner of Fig. 38, i.e., ina 
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lower power coefficient range than shown here and in the referred 


rpm range of about 350. 


6. This turbine designed for high head coefficients can now be 
compared to the theoretical values predicted in the first section of 
this thesis. The head coefficient is defined in eq. (1). The 
isentropic velocity C. is given by reference to the AT,. of Fig. 24 as 

Zz = 
C, Z /-(A)* 
| fy 


29T&p (60) 


For pressures and temperatures at the design point 


Z 
Ce . aapeeaay = 16.2 (1) 


4, UlYr 34.144 


* 
Referring now to Table I at A, = 1S: r = 0.0 and Ks = 16.0 the 
values determined from the theoretical conditions can be compared 
to the turbine design. The comparison is shown in Table IX in the 


columns Theoretical #1 and Turbine Design. 


The ratio V2 ml was arbitrarily selected equal to one for the 
theoretical development. This was not a selected parameter for the 
turbine design and the discrepancy between the two values of V2! Val 


are the basis for the differences seen in Table IX. Since the-ratio 


wA/V 


mn? was not required to be equal to one in the turbine design, 


ml 
the flow is deflected more in the rotor. This results in higher 
absolute values for es and or Since the flow deflection is 
greater, the peripheral component of the absolute velocity (Vo) is 


greater and the power output and efficiency are increased. Otherwise, 
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the theoretical predictions are good. The Me Nes value from the 
turbine design can be entered into the equations and the corres- 
pondence should be more exact. This was done and the resulting data 


is the column headed Theoretical #2 in Table IX. 


Table IX 


COMPARISON OF THEORETICAL AND TURBINE DESIGN VALUES 


Theoretical #1 Turbine Design Theoretical #2 


75" | 
0 i} 0 
1.0 
784 
955 


74.8° (after-expansion) 





* 
selected, not computed 


~61- 





7. The curves of Fig. 26 are show again in Fig. 39 with all 
values referred to the design point. All values are then shown as 
a percent of the design point values. Fig. 39 demonstrates the strong 
effect of rpm on power output and efficiency. At 70% design rpm, it 
is not possible to develop design power output. At about 90% 
design rpm, the — power can be developed at 95% of design 
efficiency. A maximum power of about 125% of the design value is 
the maximum which can be developed at 90% design rpm. At the design 
rpm, 125% design power can be developed with a loss of efficiency 
of 10%. When greater power increments are required, the efficiency 
loss becomes increasingly large. A 40% increase in power can be 


developed but at only 72.5% of design efficiency. 


At design rpm, 125% design power requires an..overall pressure 


ratio of 3.5, and 140% design power requires an overall pressure 


ratio of 8.5. 


8. As this thesis was being published, it was determined that 
the computer program finds the off-design performance of a turbine 
in which impulse conditions are maintained until choked. flow occurs 
in the rotor. A rotor blade height of 1.648 inches is required to 
satisfy such a condition. When this new rotor blade height is 
introduced into the computer program, the only direct influence is 
on the blade height correction factor for loss coefficients. The 


performance computations change is negligibly small. 
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VI. Conclusions and Recommendations 


A. Conclusions 

1. If a requirement for a high shaft power output turbine is 
encountered for which weight, size and space limitations preclude 
having several stages operating at peak efficiency, then a slightly 
higher head coefficient and slightly lower efficiency compromise 
design should not be made. Designing for power outputs higher than 
optimum should be based on extreme power output per stage. The 
initial loss of efficiency is great when a slightly higher than 
optimum head coefficient is considered. Not much more efficiency is 


lost by imposing extreme values of head coefficient. 


2. From the available literature, no straightforward method 
of computing losses is found which correlates well with the available 
data. This appears to be a fruitful area for future research and 
testing. 

a. The available data indicates the tip clearance loss is 
not zero even for impulse conditions and no clearance gap. The 
trailing vortices induced at the ends of the blades introduce 
velocity components which do no useful work. Clearance gaps within 
the annulus boundary layers should give tip clearance loss values 
about the same as zero clearance. As the clearance gap is increased, 
the loss increases. Also the losSs increases as the degree of reaction 
is increased. A high degree of reaction means a high pressure drop 
across the blade row and such a pressure differential is conducive 


to high leakage across the blade tips. 
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The loss due to tip clearance is also affected by the relative 
size of the gap to the blade height. The effect of leakage across 
the tips of the blades on the secondary flow will make measurements 
difficult to separate into tip clearance and secondary flow losses. 

b. The profile loss is a function of blade thickness to 
chord ratio, spacing, flow deflection and probably general blade 
shape. This loss offers the most possibilities for simple and 
accurate measurements in a cascade test rig. 

c. The trailing edge thickness losses are probably in- 
sensitive below a certain ratio of trailing edge thickness and 
spacing, and then assume a real importance at higher values. These 
losses, too, should be relatively simple to measure. 

d. The secondary loss measurement offers the most difficulty. 
No known theory accurately predicts the loss. Indications are that 
the secondary loss is primarily dependent on the flow turning angle 
and ratio of blade height to the chord or axial width of the blade. 
The tip clearance effect is also of primary importance and will be 
difficult to divorce in the measurements. 

e. Once the four component losses mentioned above have a 
firm basis in measurements, it will be appropriate to consider 
Reynolds number and Mach number effects as they affect the losses in 


turbine blade rows. 


3. The computer program for the turbine performance analysis 
yields results as accurate as the loss data available. The program 


only requires inputs of loss coefficients in the specified form and 


aide 





a design based on the flow function D for area determinations. 

The analysis of the off-design performance of this particular 
turbine indicates that large gains in power output can be obtained 
at pressure ratios two or three times the design value if a 15% drop 
in efficiency can be accepted. Attempting to obtain the last 
possible increment of power, however, requires extreme values of 


pressure ratio and another large loss in efficiency. 


4. Further investigation into the after-expansion behind a 
blade row is needed. A cursory examination of the magnitude of the 
axial component of the Mach number of the flow leaving an after- 
expansion showed it ot be less than one when the loss coefficient 
decreased to zero. The exact reason for this loss coefficient be- 


havior is not known. 


B. Recommendations 
1. That an organized effort be made to measure and define 

component losses in blade rows in a form useful for design purposes. 

a. It would appear that trailing edge thickness effects 
and profile losses could be measured in cascade test rigs with a 
limited number of profiles, solidities and blade heights. Some 
correlation of these measurements could be made with Ainley and 
Markov. 

b. The tip clearance loss should be measured in a rotating 
turbine test stand. The direct effects would not be hard to measure 


if the effect on secondary loss was not present. Since this effect 
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is present, it will require study to obtain quantitative correlation 
and isolate the tip clearance loss from the secondary flow effects. 
c. The secondary loss will probably have to be deduced by 
Subtracting the outer component losses from the overall efficiency 
as measured for the complete turbine blade row. No method is pres- 


ently known for measuring the secondary loss directly. 


2. Tests on the after-expansion behind a blade row should be 
made. Physical measurements may shed some light on the reason for 
the loss coefficient reaching a peak and then decreasing to negative 
values as the expansion is increased. Further theoretical study is 
also required, especially to consider three-dimensional effects. 


The testing program should include this variable. 
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FIG. / 7 Secondary losses in turbine blade rows. 
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FIG. 23. 
EFFECT OF BLADE HEIGHT 
ON LOSS COEFFICUNTS 
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FIG. 24. 
T-Ss DIAGRAM OF EXPANSION PROCESS 
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VELOCITY CT FOR TURBINE FLOW 
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DIMENSIONS OF ROTOR BLADES 
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TURBINE PERFORMANCE 
ne VALUES AND .CONSTANTS COMPUTED 


FLOW RATE = 100.0 LBM/ SEC PO = 1000. PSIA To = 1460s 2 
REFERRED FLOW RATE = GAMMA = 1.37 R = 421.50 
Cl = 326195 -€XP] = -.2701 
C2 = 318.9026 EXP2 = 3.7027 
C3 =  .000010 EXP3 = 2.7027 
Ch = 136.3089 RTO = 38.21 
PHYSICAL PARAMETERS 
RPM = 13000. ~ | REFERRED RPM = 340.226 
MEAN DIA AT ROTOR INLET = 23.00 IN STATOR BLADE HEIGHT = 1.165 IN 
MEAN DIA AT ROTOR EXIT = 23.00 IN. ROTOR BLADE HEIGHT = 1.865 IN 
STATOR DISCHARGE ANGLE = 0 DEG j 
ROTOR INLET ANGLE. _ 2 69. g9 DEG 
ROTOR DISCHARGE ANGLE’ = -71.57 DEG | 
LOSS COEFFICIENT THROUSH ROTOR (FLOWS) = .1572 
STATOR EXIT PLANE 
M(VE) VE/RTO PE/PO TESTO 
00 3 =©121.515 + 5380 +2. 8530 
STATOR DISCHARGE PLANE 
M(V1) VIJZRTO PI7PO TIVTO 3 ALFA ZETAI 
206 131.612 .4600 .8275 7h.78 .0883 
M(W1) WI/RTO PRIJPO TRIVTO  BETAI 
~80 99.071 6955 «9253 69.60 ‘ 
IN ROTOR ENTRY 
M( WI) WI/RTO PRI/PO TRIVTO BETAI 
280 O71 2.6955 .9253 69.6 
FLOW FUNCTION COMPUTATIONS 
PHI COMPUTED = .58953 PHI CRITICAL = .60837 
PRESSURE RATIO IMPOSED = 1.512 CRITICAL PRESSURE RATIO = 1.827 
ROTOR EXIT PLANE , | 3 
M( WE) WE/RTO PE2/PO TE2/TO 
of3 90.952 e 4600 e834 29 
ROTOR DISCHARGE PLANE 
‘MUW2) W2/RTO P2/PO0 T2/TO 4BETA2  ZETA2 
06 84.389 .4600 .8543 <-71.57 .27hu 
M(V2) V2/RTO P2/PO =—T2/TO AL FA2 
RESULTANT POWER AND EFFICIENCY 
VUILET/SEC) vuz FT/SEC) P2(PSIA} reetos tt T2 T2T 
. 6852.6 bekd 460 _ 3247. 1288. 


OVERALL PRESSURE RATIO = 2.17% HP = 48682. EFFICIENCY = .723 
FIG. 3 
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FIG, 32 
A-.2 DIAGRAM FOR POLYTROPIC EXPANSION - 


-97- 





Ed 


ENTHALPY 


- 


TEMPERATURE 





ENTROPY kL 


FIG. 34 
La DIAGRAM FOR AF TER-EXPANSION PROCE 55 
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Appendix I. Development of Flow Function Relations 


In Fig. 32 is shown the adiabatic expansion process in a 
turbine from state 1 to state 2. 

The fluid proceeds from state 1 to state 2' and the isentropic 
enthalpy drop is represented by process (1-2; .). The fluid proceeds 
from state 2' to state 2 and the isentropic enthalpy drop is process 
(2'-2"). For the complete process the isentropic enthalpy drop is 
represented by process 1-2... However, the sum of processes (1-2; 2) 
+ (2'-2") is greater than process (1-2...) due to the so-called 
"reheat factor''. The frictional effects raise the temperature, and 
therefore the enthalpy, of the fluid as it proceeds from state 1 to 
state 2'. This increase of enthalpy is then available to do work. 

As a result, it is necessary to make a distinction between the 
loss coefficient - for the overall process of Fig. 14 and the loss 
coefficient 4 for the polytropic process from pressure (p) to 


(p-dp) on Fig. 32. 


The polytropic loss coefficient is 


Ene] - %, (1,1) 


P 


Looking at the expansion process in more detail in Fig. 32, 


the efficiency can also be seen to be 


_ a7 
/p | 2a GE2) 
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From this the relation along the polytropic line is 


Aw 


—".~s | - (t54)* a at 0 








Ce (1,3) 
and 
—’ =<? <c, = o | I +1 
or 
AT. 7 e-/ “4. 
— Lm P (1,4) 
Since, for an adiabatic process 
B=! 
pu constant, then J c& 7 7 - (eS) 
and for an isentropic process 
* (et 
~~ = constant, then 7 LP a (1,6) 
Comparing these exponents, there is 
te (Ae) = O- ON) 
and 
== = (1,7) 


9, (671) t7 
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To find the relation between the overall loss coefficient and the 
polytropic loss coefficient, eq. (1,5) can be written between two 


points on the polytropic expansion line as 


ge 
r a a 


EG 


pe 








and 





B-T= hs & 
/s 
(1,8) 
Solving equation (1,8) for 7 and using 
p= /—4 


the reduced result is 


(Z 
4B 





To define a flow function for the mass flow through the blade 





- = e- 
Y= 





row, the expression for mass flow rate must be expressed in terms of 


the pressure ratio and the polytropic exponent. The mass flow is 


(1,10) 
“a 


=06- 





and following the procedure outlined in Appendix IV for finding the 


velocity corresponding to a pressure drop the result is 


2 ce 


yo 24 GK /—(Z = 
fy 





(IV,4) 
Along the polytropic line, for an adiabatic process 


/" 
“ee = constant 


aL 
ae ef Ne 
M5 7 qa) 


When equations (IV,4) and (1,11) are substituted into (1,10) the 


or 








result can be arranged in non-dimensional form to define a flow 


function Dp to be 


This flow function reaches a maximum at the choked condition for the 





(1,12) 


blade row. To find the corresponding pressure ratio for choked flow, 
the derivative of the flow function with respect to the pressure ratio 
will be set equal to zero and the result solved for the pressure ratio, 
2 xz / n+] 
i aS 
e- 21 -e ee 
oe / fr 
ft 7 
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Then 





2 
2 -£-/ _ tl 
aC) | £E Sy - =e ALS at 
Mr Cus] Ce ae 
which reduces to 
Ltn 
2 a {sat 
ag) _ cas yp, “Ws 
OU S-/ He Te 
When set equal to zero, the result is 
anal 
Ae = 2(_= 
re \ n#t/ 
and 
rae 
~ (2 yr 
7° CA m+] 
(1,13) 
(Tal): 


When the relation for the critical pressure ratio, eq. 
is substituted into eq. (I,12) the critical flow function is obtained 


as 


(1,14) 
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Appendix II. After-expansion from a Blade Row 


To find relations to express the expansion of a sonic or super- 
sonic flow out of a blade row, the following model is employed. The 


cascade is assumed to have infinitely many blades of zero thickness 


arranged as in Fig. 33 below. 


L 
Zz 
ON ek-04) 


Fig. 33 Relations for After-expansion from a Cascade 


A unit height perpendicular to the page is assumed. Also the uniform 
anath 
adiabatic flow is assumed to be frictionless at (w) so that R is 


perpendicular to (w) ; The control area under consideration is 


bounded by the surfaces Cw) 5 (e) and (4) » 1l.e., the "wall", the 


"exit" plane and the "discharge" plane. Ve, is sonic or supersonic. 


=109- 





The continuity equation is used to find a relation between the 


densities, 


~ [a Vv Ge (X-AX) = Ms 


4- = Aa Coz x 
and so 
SA P Ve coe (II,1) 
/e Vy] Loa («- 4X) 


The momentum equation, Ref. 1, ignoring shearing stresses is 


R =m; V, eG Vi — A, fat —r PA «+ (Te? 
where 
wi A 
we ee 
ae A . 
Ve = 2 Ye eoaan — 7 Vi den AK 
A A 
mn, = "AHA 
A AA 
my = 2 Coax — so mile el 
and so 
A 
pe =hang, - VY, (Ecoasor - fsa) tig ¢ - Ja (Etre -7' ser) 


Of particular interest are the relations in the i direction 


o = mV, -m.Yy Coz ZoK + fy © - fy a Lon x 


— a 


~*~ 





and using 


reduces to 


Le (2 
Since 
a 
a. ¢ ¢ 
*M,* |! — VY as ee Ze aay (II, 3) 
y 7Z. 


Since the process is assumed adiabatic, the total enthalpy remains 


constant through the control area, 


“x Zz 
4th 7 H+ Ke 
Z Z 
and 
Se eee ee ee 
ee ep 


and substituting 


a 
Zea +My . 2400 +V 


S—/ “2 S-/ wl. 
( p ane 


= - 





Rearranging and substituting in eq. (II,1) and using 


fa¥ = 4 fh 


results in 


2 
Ku. Me| oe tI 7 (7 Me 
Jz Vy [Low (X40) ee ( Ze) 





= (II, 5) 


Then equating eqs. (11,3) and (11,5) and collecting terms results 


in a quadratic equation in velocity ratio, 


(May 9M, cor hx — Kor X (4 ne} 


V, Loa (X-Lx) 


-(4) eo 4 / ae ee (+t) | - 
V, Crd Tee ca 6) 


The loss coefficient associated with the flow through the 


control area can be defined to be 
a 
“a — Ma 
A Z 
Vo 
based on the relations in Fig. 34. Using perfect gas relations and 


Fig. 34 as a reference, the following relations are obtained for the 


loss coefficient. 


-Ti2- 























~ 


x 


oh 


G 


= 





“iB. 
~ 
+. 
As 
| 
yen 
x 
ll 
N iS. 
~ 
~Y 
» 
| 
<Y 
pe 





(irs) 


-114- 





The Mach number of flow behind after-expansion is then, using Fig. 34 


as reference for relations, 


om, - Ve 4) ay 
AY Ye 4 |\ ALL (11,9) 
z 
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Further 
z 
i i 
Te 
az 
z& 
Mi = Me 
2 
ay 
xz 
ay. 
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= 





= ’ 


; 
® - 
* 


Substitute into (I1I,9) above 


y 
“\ a . 
ee 


i, 
Vy, we: VY, \? 
< /+ _— = (11,10) 


This computation is included in the computer program but is not 


printed out. It is inserted for ready use if a check is needed. 
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Appendix III. 


MAIN PROGRAM FOR TURBINE PERFORMANCE ANALYSIS 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 


SUBROUTINE AFTER 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 


SUBROUTINE TRNGL 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 


SUBROUTINE BEFORE 
FLOW CHART 
PROGRAM NAMES 
PROGRAM LISTING 


SUBROUTINE CURVE 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 


SUBROUTINE CPHI 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 


PROGRAMS TURPLOT AND TURBINE (HIGH HEAD COEFFICIENT DATA) 
FORTRAN NAMES 
PROGRAM TURPLOT LISTING 
PROGRAM TURBINE LISTING 


PROGRAM FANNO 3 (FLOW FUNCTION AND ZETA POLYTROPIC) 
FORTRAN NAMES 
PROGRAM LISTING 


PROGRAM FANNO 4 (CHOKED FLOW VALUES OF FANNO 3) 
FORTRAN NAMES 
PROGRAM LISTING 


PROGRAM FOR AFTER-EXPANSION COMPUTATIONS 


FORTRAN NAMES SAME AS "AFTER" ABOVE 
PROGRAM LISTING 


= bi 








- i—- 
= _ —— 
* -)' aa 


—_— ae LL 
— aq & 
-—'\|' ae 

—? = (ap 





>_> —* ee 
- '- j am 





Commute corsts 1 Gomccea | 
e: rit plane conditions 





Ae, - ee 
is < 2 = 
| Pe 
ie ~ - es ” ° 
poe isp, 7 Pw ? Subroutine AFTER 
uk caylee £ . 
{fter-expansicn co7agigagn 


i_) 
© 


els 


C44 % 5 inh | Va 
Bupro ti ee | vis fie ¥ > 0 2 


ea20ch is n ~ ya 


Suoroutine 





Subroutine TRAIGE 


Conpucc ww, , S.5 Vee 
1 al 1 








Commute TrOovor eiuay 
CON a ozs 
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PROGRAM PERFORM 
FORTRAN NAMES , EQUIVALENT SYMBOLS AND DEFINITIONS 


AE2 Ano Area at exit plane of rotor 

AJ afl Mechanical equivalent of heat 

ALFAB1 ee Undeflected flow angle of absolute velocity 
leaving stator 

ALFA1 a, Flow angle of absolute velocity leaving 
stator 

ALFA2 x, Flow angle of absolute velocity leaving 
rotor 

AME M, Mach number of flow at exit plane of stator 

AMIS M. Mach number for isentropic expansion 

AMV 1 Mo Mach number of absolute flow leaving stator 

AMV2 Myo Mach number of absolute flow leaving rotor 

AMW1 Mo Relative Mach number of flow entering rotor 

AMW2 M9 Relative Mach number of flow leaving rotor 

AMWE Mie Relative Mach number of flow at exit plane 
of rotor 

AMWEX Average Mach number of flow at exit plane 
of rotor for choked flow 

BETA1 A, Flow angle of relative flow entering rotor 

BETA2 p2 Flow angle of relative flow leaving rotor 

BETAB1 ps1 Relative flow angle for zero incidence 
entering rotor 

BETAB2 Pez a relative flow angle leaving 

GP = Specific heat at constant pressure 
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Cl 
G2 
C3 
C4 
DM1 
DM2 
DRPM 


DTE 


DTF 


DTR 


DTRAT 


DTE2 


DT2 


EFF 


EXPl 
EXP2 


EXP3 


GAM 


De 


=) 


A RPM 


(AT./T,) 
AT/T, 


(TR y7T)) 


(AT/T) 
(4T,,/T,) 


(AT,/T.) 


7 


KR 
Kus 


(Constant #1) R/g 


¥(2 & gR)/(4 -1) 
(@-1)/(2 & gr) 


Y Mor 


Mean diameter of blade at rotor entrance 


(Constant #2) = 


(Constant #3) 


(Constant #4) 


Mean diameter of blade at rotor exit 
Increment of rpm 


Temperature drop from inlet to exit plane 
of stator 


Final temperature difference from total to 
static temperature after rotor 


Temperature drop between relative total 
temperature ahead of rotor and static 
temperature ahead of rotor 


Temperature drop from inlet to discharge 
plane of stator 


Temperature drop from inlet to exit plane 
of rotor 


Temperature drop from inlet to discharge 
plane of rotor 


Efficiency based on total pressure after 
rotor 


Exponent #1 = (®@-1)/*% 
Exponent #2 = %/( %-1) 
Exponent #3 = 1/( ®-1) 


Acceleration of gravity 
Specific heat ratio 


Blade height correction for rotor loss 
coefficient 


Blade height correction for stator loss 
coefficient 


Sea 





HP 


Hs 


PE 


PE2 


PRC 


PHICR 


PPCR 


PO 


PPRAT 


PRAT 


PRATO 


PR1 


PRIY 


Pl 


PZ 


P2T 


RFR 


RPM 


S1 


(p,/P) 
(P.o/P) 
p 

Pex 


Gi®) + 


Pei /P2 


P,/P, 


Po/Po 
Fai’? 
mao 
aa 


p,/P 


p/P 


W TO/P 


RPM/ } Lid 


Horsepower 
Rotor blade height 
Stator blade height 


Number of data points stored for a 
particular loss coefficient curve 


Static pressure at exit plane of stator 
Static pressure at exit plane of rotor 
Flow function computed with input % and ss 


Critical flow function for choked flow in 
rotor 


Critical pressure ratio for choked flow in 
rotor 


Total pressure ahead of stator 
Pressure ratio through rotor 


Pressure ratio between stator exit plane 
and stator discharge plane 


Overall pressure ratio 
Relative pressure ratio 
Relative total pressure ahead of rotor 


Relative total pressure behind shock in 
rotor entry 


Static pressure between stator and rotor 
Static pressure after rotor 

Static pressure after rotor 

Gas constant 

Referred flow rate 


Referred rpm 
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RPMA 


RTO 


TEIS 


TI 


TO 


Tl 


1074 


21 


TP2 


TRAT 


TRATR 


TR1 


rrZ 


m1Y 


Ul 


U2 


UV 


vul 


Tae) T 
eo 
Cis 


e2/T 


Actual value of rpm 


Temperature at exit plane of stator 


Isentropic temperature at exit plane of 
Stator 


Temperature at exit plane of rotor 


True interval between data points on loss 
coefficient curves 


Total temperature ahead of stator 

Static temperature between stator and rotor 
Static temperature after rotor 

Static temperature after rotor 

Total pressure after rotor 

Isentropic temperature ratio through stator 
Isentropic temperature ratio through rotor 
Relative total temperature ahead of rotor 
Total temperature after rotor 


Static temperature behind shock in rotor 
entry 


Peripheral speed of blade at mean diameter 
at rotor entrance 


Peripheral speed of blade at mean diameter 
at rotor discharge 


Abscissa for rotor loss coefficient data 


Absolute velocity of flow at exit plane of 
Stator 


Meridional component of absolute velocity 


Peripheral component of absolute velocity 
leaving stator 
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vu2 


V1 


V2 


Wl 


W1Y 


W2 


ZET 


ZETAA 


ZETAB 


ZETAS 


ZETAA1 


ZETAA2 


ZETA1 


ZETA2 


Vu2/ die 


V,/% 2, 
v,/ 7 um 
Ww 


a. ] To 
Ww! y - 
Waa! Y E 


Wao! Yt 


O 
Wt! " 


oO 


Peripheral component of absolute velocity 
leaving rotor 


Absolute flow velocity leaving stator 
Absolute flow velocity leaving rotor 

Mass flow rate 

Relative velocity at exit plane of rotor 
Meridional component of relative velocity 


Peripheral component of relative velocity 
entering rotor 


Peripheral component of relative velocity 
leaving rotor 


Relative flow velocity ahead of rotor 


Relative flow velocity behind entry shock 
in rotor 


Relative flow velocity behind rotor 
Ordinate of loss coefficient curves 


Loss coefficient for velocity through rotor 
including after-expansion 


After-expansion loss coefficient 

Loss coefficient due to rotor entry shock 
Stator loss coefficient 

Loss coefficient for flows through stator 
Loss coefficient for flows through rotor 


Loss coefficient for velocities through 
stator 


Loss coefficient for velocities through 
rotor without after-expansion 
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ALFA 


COSRAT 


DALFA 


DALFAR 


PRAT 


PRATI 


VD 


VELRAT 


ZA 


SUBROUTINE AFTER 


FORTRAN NAMES, EQUIVALENT SYMBOLS AND DEFINITIONS 


O*;, 


Ac&K 

AX 
(p4/P.,) 
(Py/P.)4 
V 


d 


Ae 


ys 


V 


Undeflected flow angdl (blade angle) 
Cosine ratio of blade angle to flow angle 
oO, ~- X. (degrees) 
B lL 
of, x, (radius) 
Imposed pressure ratio 
Computed pressure ratio 
Velocity at discharge plane 


Velocity ratio of flow at discharge plane 
to exit plane 


Loss coefficient due to after-expansion 
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AbsolLuteniuiciwes 
wind, rerrolwe 
) = 2 BPOLotive 421 eer ou 
j sbpsolute 





SUBROUTINE TRNGL 


FORTRAN NAMES, EQUIVALENT SYMBOLS AND DEFINITIONS 


A Input angle of flow 

B Output flow angle (If A= ax ; B= 8 and ifa= 2 ; 
B=.) 

N Value selects plus or minus SIGN 

U Peripheral velocity of blade at mean diameter 

VM ve Meridional component of absolute velocity 

VU1 Vil Peripheral component of absolute velocity 

WUl Wa Peripheral component of relative velocity 

Wl W, Relative velocity 
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ALFA 


ALFAI 


CC 


DALFA 


DTR 


FF 


GG 


PP 


PRP 


PR1Y 


Pay 


ZETAB 


_SUBROUTINE BEFORE 


FORTRAN NAMES, EQUIVALENT SYMBOLS AND DEFINITIONS 


AX 


(T4°T)) 


P,/P,* 
Pai/Py 
faa Fo 
Bie 


os 
Ws MM, 


Interim computation for velocity ratio 
Blade angle entering rotor 
Flow angle ahead of rotor 
Mach number of flow after shock 
Interim computation for velocity ratio 
Cosine ratio of angles involved 
2 ‘ 
A - B, used to determine whether shock should 


be computed, or the component of the velocity 
in direction of blade be used. 


AK = XK, - x, 


Temperature difference between relative total 
temperature and static temperature 


Interim computation to find Mach number after 
shock 


Interim computation to find Mach number after 
shock 


Static pressure ratio through shock 
Pressure ratio through shock 
Relative total pressure behind shock 
Static pressure behind shock 


Velocity ratio of flow ahead of and behind 
shock 


Loss coefficient of flow due to entry shock 
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~~. .., = Nis cts 


' 
Boe 
i” i wows 















, OURPOT 
SOmege Ord mere 


Celi ih array (3) 





Write equation of paravola 
throuch three nearest points 


Compute ordinate 


moo. eet 











— 


N 
NMAX 


ORD 


gh 


T2 T1/2 


FORTRAN NAMES, EQUIVALENT SYMBOLS AND DEFINITIONS 


SURROUTINE CURVE 


> 
The terms of equation ORD = a + bx + cx 


Subscripts tor ordinate data points 


Counter to locate three nearest data points 
to input abscissa 


Array selector 
Number of data points in selected array 


Ordinate computed, i.e., loss coefficient 
value desired 


Data abscissa interval 


One-half of data interval. Used to select the 
three nearest points 


Abscissa in terms of three nearest points 


Input abscissa 
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CPHI 
DAB 
DAN 


DIFF 


EN 

EXP1 
EXP4 
EXP6 
EXP/ 


GAM 


PHI 2 


PHICR 


PPCR 


RC 


RB 


ZETA 


ZETAP 


KX & 


p 


cr 


(P/p) 


Sp 


SUBROUTINE CPHI 


FORTRAN NAMES , EQUIVALENT SYMBOLS AND DEFINITIONS 


Interim critical pressure ratio while computing 
Compute D , Subroutine name 


Name for the term Gina! 


Name for the term Gir: 


Difference between given and computed loss 
coefficient 


Polytropic exponent 
Exponent #1 = (@ - 1)/ ® 
Exponent #4 = (1 - %) (S- Des 


2/n 


Exponent #6 


Exponent #7 = (nt1)/n 

Specific heat ratio 

Flow function computed 

Flow function, squared 

Choked flow value of flow function 
Choked flow value of pressure ratio 
Computed loss coefficient 


Name for the term (pee 


Name for the term (peyote 


Loss coefficient 


Polytropic loss coefficient 
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AKE 
AKIS 


ALFA1 


ALFA2 


BETAL 


BETA2 


DBETA 
DETA 
EFF 
ETA 
ETAA 


ETAC 


EXX 


NALFAI 


NPHIE 


PHI 


PHIR 


PSI 


PROGRAMS TURPLOT AND TURBINE 


FORTRAN NAMES , EQUIVALENT SYMBOLS AND DEFINITIONS 


% 


Rie 


«< 6x BF 2” 


Leaving loss coefficient 
Head coefficient, eq. 1 


Angle of flow exiting stator (radians), 
absolute 


Angle of flow leaving rotor (radians), 
absolute 


Flow angle entering rotor (degrees), 
relative 


Flow angle leaving rotor (degrees), 
relative 


a6 = 2, - P2 

a 75 

Reheat factor, defined in Ref. l 
Efficiency 

Efficiency assumed 


Efficiency computed as an interim step for 
comparison 


Expression defined by eq. 5a 


Angle of flow exiting rotor (degrees), 
absolute 


Carryover coefficient of kinetic energy into 
next stage 


Velocity coefficient in stator 


Carryover coefficient of kinetic energy from 
Stator into rotor 


Velocity coefficient in rotor 
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‘ ——_—)' & & &= 4h ea ae . 
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PSIA 


RADRAT 
RBETA1 
RBETA2 
RELRAT 


RSTAR 


VELM1 


VE LRAT 
WSQR1 
WSQR2 


XSQR 


Ro/R, 


Velocity coefficient in rotor computed for 
comparison 


Radius ratio through rotor 

Flow angle entering rotor (radians), relative 
Flow angle leaving rotor (radians), relative 
Ratio of relative velocities 

Degree of reaction, that fraction of the 
expansion through the stage which occurs in 
the rotor 

Meridional component of absolute velocity 
leaving stator/Peripheral velocity of rotor 
at mean diameter 

Velocity ratio of meridional components 

The relative velocity entering rotor, squared 


The relative velocity leaving rotor, squared 


Expression defined by eq. 5a 
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DAB 
DAN 
EE 
EN 
EXP1 
EXP2 
EXP3 
EXP4 
GAM 
PH1 


PH12 


RB 


ZETA 
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PROGRAM FANNO3 


FORTRAN NAMES , EQUIVALENT SYMBOLS AND DEFINITIONS 


P/p 


6 B® 3 


Pressure ratio 


Name for the term (Gey 


Name for the term qin 


Reheat factor defined by ( - ¥y/a 5): 
Polytropic exponent. See Appendix I 

Exponent #1 = (2/n) 

Exponent #2 = (ntl)/n 

Exponent #3 = (1 - ¥,) ( @- 1)/® 

Exponent #4 = (4% - 1)/% 

Specific heat ratio 

Flow function 

Flow function squared. Intermediate step in finding 


Name for the term (p/p) e*P3 


Name for the tem (aneare 


Computed loss coefficient to be used for 
determining p Rk (RA-RB)/(1-RB) 


Loss coefficient # 1 -— 


Polytropic loss coefficient defined in Appendix I 
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PROGRAM FANNO4 


FORTRAN NAMES , EQUIVALENT SYMBOLS AND DEFINITIONS 


(P/P). 


p 
gr 
B, 


Pe. 


(P/p) 


P/p 


cE 


rc 


Interim critical pressure ratio while computing 
Compute D » Subroutine name 

Name for the term (piper? 
Name for the term (p/pye*P> 


Difference between given and computed loss 
coefficient 


Polytropic exponent 

Exponent #1 = (& - 1)/@ 

Exponent #4 # (1 - ¥») (os - 1)/% 
Exponent #5 = 2/n 

Exponent #6 * (nt1)/n 

Specific heat ratio 


Critical pressure ratio name for storing away 
and printing 


Flow function assumed for comparison 
Flow function, squared 


Flow function computed in subroutine based on 
® and 


Choked flow value of flow function 
Choked flow value of pressure ratio 
Pressure ratio 

Computed loss coefficient 

Name for the term (p/p) °*P4 


Name for the term (p/p) Pt 
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Appendix IV. Critical Mach Number of Flow Undergoing Polytropic 
Expansion. 


The critical pressure ratio for a flow undergoing a polytropic 


expansion, as developed in Appendix I, is 


nae. 
—) - (4 
? ZOn8 & (1593) 


To find the maximum velocity corresponding to this pressure 


ratio it must be recalled that, for a polytropic process 


ww 
7 ww = constant 


(IV,1) 
Tr s-7t m 
Pre = r y = constant 
na! 
Tr 
7=(C 7 PS 


where C is a constant 


Therefore 











r= ofa) 
- (IV, 3) 


Using the T-s diagram of Fig. 24 as a reference, it can be seen that 


m= 
,-T = 7 |/ tH a 
, 


~lio7- 


hea 


rot 





and 








zZ a 
a. = O(%-7) = £*t % |/ -(£) - 
hgJ J Y-! Fe 
(IV,4) 
For the critical velocity then, using (1,13) 
Zi n=l 
Yr = 27 * BI/- / ” 
ice! (2 = 
x 
which reduces to 
Y~ _ 2gkT ¢ (21) 
7 See 
T-1 \ mei 
(IV, 5) 


The acoustic velocity is 


Z ad Ye!) 
a’ = 69kT = per (Z) . a7 (at 
p Z 


which reduces to 








me” = wget ( Z ; 
Therefore 
Zz 
M2 Ma 2 el 
Zz 
ta a 7 ante 
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Also developed in Appendix I was the relation for the polytropic 


exponent n, in terms of the polytropic loss coefficient, as 


oO 
Xp (4-1) +/ 


fe = 





(1,7) 


Making this substitution, and collecting terms 
’ 
fet sey / —%p 
/ Y —-S-/ 
ts % (IV,8) 


and 





Ma ie 
i /+ 340 -/) 


(IV,9) 


For example, for Oo = 1.37 and ¥ = 0.10, M = .931, 
p er 

and for ¥ = 0.20, M_ # .854. 
p er 
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